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S.  M.  Lee  and  C.  B.  Bell 
San  Diego  State  University 


1.  Introduction  and  Summary. 


NHPP's  (i.e.,  non-homogeneous  Poisson  processes)  describe  a  wide 
variety  of  physical  phenomena-  background  noise,  epileptic  seizure 
patterns,  failure  patterns  etc.  A  particularly  useful  NHPP  is  the 
Weibull-type  NHPP  -  one  for  which  the  mean  function  is  of  the  form 

g 

y(t)  =  at.  It  is  immediate  that  by  appropriate  choice  of  B,  the 
shape  parameter,  one  can  achieve  a  constant  failure  rate,  and  an 
increasing  failure  rate,  as  well  as  a  decreasing  failure  rate.  Con¬ 
sequently,  this  family,  ft(WPP),  of  Weibull-type  NHPP's  incorporates 
a  wide  variety  of  useful  models. 

The  developments  in  the  sequel  will  be  concerned  with  interval 
and  point  estimation,  as  well  as  signal  detection,  which  is,  of  course, 
closely  related  to  hypothesis  testing. 

In  deriving  the  inference  procedures,  one  will  consider  four 
different  sampling  plans,  which  will  be  described  in  detail  in  the 
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sequel.  They  are  (A)  Type  I  Censoring;  (B)  Type  II  Censoring, 

(C)  Same-shape  Sampling;  and  (D)  Same-distance  Sampling. 

It  will  be  found  that  the  requirements,  the  data,  and  the  feasible 
types  of  inference  are  different  for  these  different  plans. 

Hie  paper  is  divided  into  twelve  sections.  Soiu  basic  definitions 
and  background  material  is  given  in  Section  2.  Section  3  contains  a 
description  of  the  sampling  plans .  Sections  4  and  5  present  the 
necessary  distribution  theory. 

Detection  problems  are  treated  in  Sections  6  through  9.  Estimation 
is  treated  in  Sections  10  and  11.  Open  problems  and  conclusions  are 
presented  in  the  final  section.  Finally,  numerical  illustrations  of 
the  methodology  developed  are  presented  in  the  appendix. 

2.  Some  Basic  Definitions  and  Distributions. 

Hie  stochastic  law,  L,  of  a  >HPP  is  completely  determined  by 
the  mean  function  of  the  process.  Mean  functions  here  will  be  of 
slightly  restricted  type. 

Definition  2.1.  A  mean  function  v(’)  is  a  real  valued  function 
(defined  for  non-negative  time  and)  satisfying  (i)  y(0)  =  0; 

(ii)  y(®)  -  »;  (iii)  y(*)  is  continuous;  and  (iv)  y C * )  is 
strictly  increasing. 

Some  common  mean  functions  are  given  below. 
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Example  2.1.  (i)  For  HPP's,  i.e.,  homogeneous  Poisson  processes. 


one  has 

y(t)  =  At, 

for  some 

A  >  0. 

(ii) 

For  each  law 

L  in 

Q(WPP) , 

p(t)  =  at6, 

for  some 

positive  a 

and  8.  One 

notes  that 

when  S 

=  1,  one  has 

a  HPP. 

(iii) 

yft)  = 

a  ln(l  +  bt) . 

(iv)  p(t) 

-ln[$(-pln  At)],  where  $(•)  is  the  standard  Gaussian  cpf. 

(v)  p(t)  =  a(eP  -  1),  related  to  the  Gumbel-Extreme-Value  density 

(vi)  p(t)  =  At  exp  (Z*  B),  where  Z  and  B  are  vectors.  This 

'Xi  'Xj  'X,  'X, 

corresponds  to  the  Cox  proportional  hazard  model 
A  NHPP  can  now  be  Formally  defined. 

Definition  2.2.  (a)  Consider  a  counting  process  (N(t):  t  O) 

satisfying  (i)  N CO )  =  0,  (ii)  increments  are  independent,  and 

Ciii)  P{N(t)  -  N(s)  =  k)  =  [v(t)  -  p(s)]k  exD  {~[p(t)  -  p(s)])/k! 

0  _<  s  <_  t ,  and  k  »  0 ,  1 ,  2 ,  ...  where  p  ( • )  is  a  mean 

function.  Such  a  point  process  is  the  counting  process  of  a  NHPP  with 

mean  function  u(Q.  (b)  {Wr:  r  ^  1}  is  the  (dual)  waiting-time 

process  of  the  NHPP  with  mean  function  u(*)»  if  the  Ws  satisfy 

(N(t)  <  n)  =  {W  .  >  t)  for  n  =  0,  1,  2,  ...  .  (c)  (X  :  r  >  l) 

—  n+i  t  — 

is  the  (dual)  interarrival-time  process  if  X  =  W  -  W  where 


There  are  some  special  relations  between  and  p(*)  which 

characterize  the  stochastic  law,  L,  of  a  NHPP. 

Theorem  2.1.  Let  F(t)  »  P{W^  £  t)  and  f(t)  =  .  Then 

(1)  F(0)  =  0  and  F(t)  =  1  -  exp  {-y(t)},  for  t  ^  0. 


'*  . 
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(2)  f(t)  =  y'(t)  exp  (-y(t)},  for  t  ^  0. 

C3)  v'(t)  =  fSctT  and 

(4)  y(t)  *  -In [1  -  F (t) ] . 

In  the  sequel,  one  will  be  interested  solely  in  ft(WPP),  the 
family  of  Weibull-type  NHPP  laws  L. 

Definition  2.3.  A  NHPP  has  a  law  L  in  ft(WPP)  if  its  mean 

a 

function  is  of  the  form  y(t)  *  atp  (t  0),  for  positive  a  and  3- 
The  other  relevant  functions  for  ft (WPP)  are  as  follows. 

Q 

Theorem  2.2,  For  L  in  ft(WPP)  with  y(t)  =  at  ,  one  has 

(i)  y'(t)  *  aBt0'1,  t>0, 

(ii)  F(t)  =  1  -  exp  {-at  },  t  _>  0;  and  f(t)  *  aBt^"1  exp  {-at**},  t  ^  0. 

The  first  waiting  time  Wj,  then,  has  a  (two-parameter)  Weibull 
distribution.  This  accounts  for  the  name  of  this  particular  family 
of  NHPP’s. 

By  an  appropriate  change  of  scale  and/or  an  appropriate  choice  of 
parameters  a  WPP  can  be  transformed  to  a  HPP,  i.e.,  homogeneous 
Poisson  process. 

Theorem  2.3.  Let  (N(t):  t  ^  0}  have  a  law  L  in  ft(WPP), 

8  1/8 
i.e.,  y(t)  *  at  for  some  positive  a  and  3;  let  M(s)  ■*  N(s  '  ), 

3 

for  s  >  0;  let  V  *  W_;  and  Z_  ■  V  -  V  . , 

—  T  T  T  T  r-1 


Then 


where  VQ  ■  WQ  ■  0. 
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(1)  (N(t):  t  >_  0}  is  a  HPP  iff  $  =  1; 

(2)  {M(t) :  t  0}  is  a  HPP  with  y(t)  =  at; 

(3)  (Vr)  and  { }  are,  respectively,  the  waiting-times  and 

interarrival  times  of  (M(t)};  and 

(4)  {Z^}  are  i.i.d.  Exp(a). 

These  insights  will  be  used  in  constructing  the  detection  procedures 
in  the  sequel. 

Also,  one  will  employ  the  Kolmogorov- Smirnov  statistic  (Kolmogorov, 

1933)  and  its  modifications  based  on  the  works  of  I.illiefors  (1967,  1969); 

Srinivasan  (1970)  and  Choi  (1980). 

Let  X  =  (X. ,  ....  X  )  be  a  random  sample  with  common  continuous 
'v,  i  n 

cpf  F(-  ;0) ;  and  let  ft"  «  {F(-;8):  0  E@).  Let  ft"  admit  a  M-S-S 

(minimal  sufficient  statistic),  S(X).  Further,  let  X*  =  [X(l),  ..., 

'Xj  % 

X(n)]  be  the  order  statistics  of  X  and  F  (*)  be  the  empirical 

%  n  "  l  n 

cdf  of  the  X's ,  i.e.,  F  (z)  =  -  £  e  (z  -  X.)  =  -  £  e  ( z  -  X(j)). 

n  n  j=l  J  1 

The  following  statistics  will  be  employed  in  the  sequel. 

Definition  2.4.  Let  0  be  an  arbitrary  fixed  element  of 

A 

and  0  be  the  MLE  (maximum  likelihood  estimate)  based  on  X, 

*\j 

(i)  Pn(z)  =  F(z;9),  for  all  z. 

(ii)  F  (z)  *  E(F  (z)|  S(X)>; 
n  n  'v 

(iii)  DnC9)=  sup  |Fn(z)  -  F(z;0)|; 
z 

(iv)  Dn  *  sup  | Fn  Cz)  -  Fn(z)|;  and 
z 

(v)  Dn  =  sup  |Fn(z)  -  Fn(z)|. 
z 
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The  cpfs  of  these  last  three  statistics  are  utilized  in  treating 
signal  detection.  The  notation  to  be  used  for  these  and  other  relevant 
distributions  is  given  in  Definition  2.5.  below. 

Definition  2.5.  Tbe  cpfs  K-S(n);  ft"  -  LI(n);  ft"  -  SR(n); 
G-O-S(n);  U-O-S(n)  are  defined  by  the  following  distributions  of 
statistics. 

(1)  Dn(9)  -x-  K-S(n); 

(2)  Dn  ^  12"  -  LI(n); 

(3)  Dn  -v  ft"  -  SR(n); 

(4)  X* -v.  G-O-S(n),  when  G(-)  ■  F(- ,0rt) ;  and 

(5)  X*  'v  U-O-S(n),  when  G(')*U(0,1). 

Some  standard  distributions  are  denoted  using  the  notation  below. 

NOTATION ;  (1)  refers  to  the  chi-square  distribution  with  m 

degrees  of  freedom. 

(2)  F(m,r)  refers  to  the  classical  F-distribution  with  m  and  r 
degrees  of  freedom. 

(3)  r(m, A)  denotes  the  gamma  distribution  with  parameters  m  and  A. 

(4)  H(-jY>A)  is  the  cpf  satisfying  H(z;y,A)  <*  (| for  0  <  z  <  A, 
with  y  and  A  >  0. 

(5)  NB(m,p)  denotes  the  negative  binomial  cpf  such  that  P(X  *  r}  * 
(m+y_1)prqm  for  r  *  0,  1,  2,  ...  . 

One  next  presents  the  four  sampling  plans. 


3.  The  Four  Sampling  Plans. 


For  life-testing  and  survival  analysis  both  in  engineering  fields 
and  in  biomedical  testing,  censoring  arises  naturally. 

Plan  I.  Type  I  Censoring. 

One  observes  the  WPP  on  the  interval  [0,T*],  and  records 
[N(T*J;  Wj,  ....  WN(T#)]  =  [k*.  Wj,  ....  »k],  where  N(T*)  -  k  >  1. 

This  sampling  plan  is  used  when  the  time  interval  available  for 
study  is  limited.  It,  of  course,  could  happen  that  very  few  events 
occur  in  [0,  T*],  i.e.,  k  is  small.  If  such  an  occurrence  has  a 
reasonably  probability,  one  might  wish  to  consider  an  alternate 
sampling  plan. 

Plan  II.  Type  II  Censoring. 

Here  one  observes  the  WPP  until  the  kth  waiting  time,  W^; 
and  records  W  =  [W^,  ....  W^]. 

This  sampling  plan  is  used  when  a  certain  minimum  amount  of  data 
is  needed.  It  could  happen  that  W^  is  much  larger  than  anticipated, 
and  one  cannot  afford  the  time.  Consequently  one  might  wish  to  consider 
some  combination  of  Plans  I  and  II.  (Such  a  combination  will  not  be 
considered  in  the  sequel). 

Plan  III.  Same-Shape  Sampling. 

This  is  a  direct  adaptation  of  a  sampling  plan  used  by  Basawa  and 


Rao  (1980)  for  HPP’s. 
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Let  {N(t) :  t  _>  0)  be  the  WPP  (of  interest)  with  mean  function 
0 

u(t)  =  at  and  waiting  times  {Wn>;  and  let  {N*(t):  t  _>  0}  be 

g 

an  independent  WPP  with  mean  function  U*(t)  =  t  ;  and  waiting 
times  {W*}.  These  two  mean  functions  have  the  same  shape,  i.e., 

[y(t)]  [v>*(t)]  1  is  constant.  This  plan  can  be  implemented  only  when 
3  is  known . 

One  observes  the  process  of  interest  only  at  the  waiting  times 

{W*}.  The  data  is  then 
n 

N  *  [N(wp,  ....  N(W*) ]  or  Y*  «=  [Yf,  Y*],  where 

Y*  =  N(W*)  -  N(W*  _)  with  W*  *  0. 
r  r  r-i  u 

Sampling  in  this  fashion  yields  tractable  distributions  for  some 
statistics  of  interest.  In  fact,  one  can  prove,  analagous  to  Basawa 
and  Rao  (1980) 

Theorem  3.1.  (1)  YJ,  ...,  Y£  are  i.i.d.  MB(1,  and 

(2)  N(W*)  =  l  Y.  %  N  -  B(k,  y^). 

Plan  IV.  Same-Distance  Sampling. 

One  observes  the  process  (N(t)}  at  times  0  <  tj  <  <  ...  <  t^, 

WO 

where  y(t  )  =  ry(t,).  This  means  that  t  =  t,(r  '  ),  for  r  *  1,  2, 
r  l  r  l 

k;  and  that  one  must  know  8,  to  implement  the  plan. 

The  data  here  is  N  =  [N(t,),  ...,  N(t.)J  or  Y  *  [Y.,  ....  Y.] 

*\j  A  K  1  K 

where  Y  *  N(t  )  -  N(t  .)  with  t_  »  0.  The  Y's,  in  this  case, 
r  r  r-i  u 


-  9  - 

constitute  a  random  sample,  and,  hence,  some  of  the  statistics  to  be 
used  will  have  more  tractable  distributions. 

The  distribution  theory  of  interest  is  given  below. 

4.  Distribution  Theory  when  the  Shape  Parameter,  6,  is  known 

When  3  is  known,  the  same-shape  and  same-distance  samplin  .s, 

which  will  be  described  in  the  sequel,  are  best  handled  by  transforming 
the  WPP  to  a  HPP  (as  in  Theorem  2.3),  with  mean  function  y(t)  =  at. 
Inference,  then  concerns  only  one  parameter  (and  only  HPP's). 

For  Type  I  and  Type  II  censoring,  one  might  wish  to  bring  in  some 
different  types  of  results. 

(A)  Type  I  Censoring. 

Here  one  observes  the  process  on  the  time  interval  [0,T*],  and 
receives  data 

[N(T*) ;  Wr  W2,  ...,  WN(T^J  =  [k;  WJf  ....  Wk]  when 
N(T*)  =  k  _>  1. 

For  such  data  one  has  the  following  useful  distribution  theorem 
for  12  (WPP) . 

Theorem  4.1.  Conditionally,  given  N(T*)  =  k(>  1), 

(i)  [W  ,  ...,  Wk]  'vG-O-S(k),  where  G(-)  =  H(-,$,T*)  (See  Section  3). 
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we  w  3 

(ii)  I tp)  .  •••.  Cp  ]  ^U-O-S(k);  and 

(iii)  [Yj,  Y2,  ....  Yk]  %  Exp (6) -O-S(k) ,  where 
W  i  • 

Y.  =  In  (■  - ftr1)  for  1  <  j  <  k. 

These  order-statistic  distributions  lead  naturally  to  the  K-S 
(and  other  goodness-of-fit)  statistics. 

Theorem  4.2.  Conditionally,  given  N(T*)  =  k(>  1),  each  of  the 

following  statistics  is  distributed  K-S(k). 

(i)  sup  |i  l  e  (z  -  W  )  -  (=y)0|, 
z  K  1  J  1 

k  W 

(ii)  sup  |£-  l  e  (u  -  fc=i]6)  -  «! 

0<u<l  K  1  1 

(iii)  sup  l  e  (z  -  Y.)  -  [1  -  e‘zg]| 
z  K  1  3 

For  Type  II  censoring,  the  two  theorems  above  change  only  slightly. 

(B)  Type  II  Censoring 

Here  one  observes  the  WPP  until  the  ktn  waiting  time  W^;  and 
the  data  is  W  =  [W^ ,  . . . ,  W^] . 

Theorem  4.3.  Conditionally,  given  ■  w*, 

(i)  [Wj,  ....  Wk  j]  n,  G-0-S(k  -  1),  where  G(-)  =  H(*,g,w*); 

W  3  W  3 

(ii)  [(57)  .....  ]  *  U-0-S(k  -  1);  and 
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(iii)  [Yj,  Y^]  Exp(S)-0-S(k  -  1),  where 

w 

Yj  =  -In  for  1  <  j  <  k  -  1. 

The  relevant  K-S  distributions  are,  then,  as  given  below. 

Theorem  4,4.  Conditionally,  given  =  w*,  each  of  the  following 
statistics  is  distributed  K-S  (k  -  1). 

1  k_1  ft 

(i)  sup  |  ~  l  c  (z  -  W.)  -  (^)P|; 

k- 1  W . 

(ii)  sup  |  l  e  (u  -  [-i]3)  -  u| ;  and 

0<u<i  1 

1  _Q 

(iii)  sup  |  t-j  l  e  (z  -  Y.)  -  [1  -  e  P) | . 

z  K_l  1  3 

If  both  a  and  8  are  known,  one  can  utilize 
Theorem  4.S.  2y(W^)  =  2aW^  'v 

When  $  is  unknown,  the  problems  are  significantly  different.  As 
a  nuisance  parameter,  the  shape  parameter  3  cannot  be  eliminated  in 
the  same  way  one  would  eliminate  the  scale  parameter  a. 

5.  Distribution  Theory  when  £  is  unknown 

When  the  shape  parameter  is  unknown,  one  can  still  construct  a  HPP 
from  the  data,  but  the  development  is  more  involved. 


12  - 


One  starts  with  Type  I  Censoring  and  Theorem  4.1,  in  which 
Y  =  [Yr  ....  Yk]  *  Exp  (B)-O-S(k) ,  where 
W 

Y.  =  -In  f  K*^].  Now  let 

r  k  w 

z*  -  I  (k  +  1  -  j)Y  =  -  l  m  ln(=£)  and 
j=l  -  m=k+l-r  1 

1* 

E  «  {Ej,  E2,  ....  EJc_1)  ■  5  1  1  r  <  k  -  l}. 

One  can  then  prove 

Theorem  5.1.  Conditionally,  given,  NTT*)  *  k, 

(1)  Z*  ^  r(r,B)  and  2BZ*  U?  for  r  =  1,  2,  ....  k; 

r  r  zt 

(2)  E-v-U-O-SCk  -  1);  and 
% 

(3)  Z£  and  E  are  independent. 

Since  neither  E  nor  its  distribution  involves  the  (nuisance) 
parameters  a  and  0,  it  can  be  used  to  investigate  questions  of  the 
structure  of  the  process . 

For  the  two-sample  situation,  one  "eliminates"  8  in  an  analagous 
fashion. 

Consider  two  processes.  The  first  {N(t)>  has  law  Lj  and  waiting 
times  W  *  (W, ,  ...»  W  )  in  the  interval  [0,  T*];  and  the  second 

^  l  m 

process  (N*(t)}  with  law  L 2>  has  in  the  same  time  interval  waiting 

time  V*  »  (V*,  ...,  V*).  Now  let 
i  n 
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W 

Yj  =  -In  (~~y  for  1  _<  j  <_  n,  and 
V* 

Uj  =  -In  C— jt'1)  for  1  <  j  <  n. 

Then,  one  readily  proves 

Theorem  5.2.  Conditionally,  given  Nj(T*)  =  m  and  NjfT*)  =  n, 

— 'v  F(2m,  2n),  if  L.  and  L-  are  in  fl(WPP)  and  B.  =  B_. 

U  1  1  1 

One  notes  here  that  a  does  not  explicitly  appear  in  the  results 
above.  However,  the  value  of  a  does  influence  the  distributions  of 
{N(t)}  and  (N*(t)} 

For  Type  II  Censoring,  the  results  are  quite  similar. 

Theorem  5.3.  (i)  Conditionally,  given  W^+1  =  w*,  E  and 

have  the  distributions  in  Theorem  5.1  if  w*  replaces  T*  in  the 
definition  of  Y^. 

(ii)  Conditionally,  given  Wffl+1  =  w*  and  Vn+j  =  v*,  the  con¬ 
clusion  of  Theorem  5.2  is  valid,  when 

Yj  =  -m  (%£=*)  and  u.  =  -m 

One  can  now  treat  the  signal  detection  problems. 

6.  Signal  Detection:  Goodness -of- fit  with  Two  Nuisance  Parameters. 
The  signal  detection  problem  here  is  as  given  below. 

PNj :  L  e  fl(WPP)  vs.  N  +  Sj  L  ^  ftfWPP) 
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Based  on  Theorems  4.1  and  5.1  and  Definitions  2.4  and  2.5,  one 
can  construct  the  following  decision  rules. 

Decision  Rule  6.1.  For  Sampling  Plan  I,  i.e..  Type  I  Censoring 
decide  N  S^  iff 

.  k-1 

Dk-1  =  sJjp  ^  e  -  U|  >  d(a,k-l) ,  where 

d(a,  k-1)  is  the  appropriate  percentile  of  the  K-S(k  -  1),  distribution, 
and  the  E's  are  as  in  Section  5. 

Decision  Rule  6.2.  For  Type  I  Censoring,  decide  N  ♦  iff 

k 

^  »  sup  l  e  (z  -  Y.)  -  [1  -  exp  {-  ^)]|  >  d(a,k) 
z  K  1  3  Y 

/s 

where  d(a,k)  is  the  appropriate  percentile  of  the  ft"- LI 
distribution  of  Definition  2.5  and  Lilliefors  (1969),  with 
ft"  =  {Exp(B):  B  >  0}. 

Decision  Rule  6.3.  For  Type  I  Censoring,  decide  N  ♦  iff 
K 

D.  -  sup  |i  l  z  (z  -  Y  )  -  [1  -  (1  -  -r-)1"1]!  >  d(a,lc) 
k  x  i  ;  ky 

'V, 

where  d(o,k)  is  the  appropriate  percentile  of  the  ft"-SR(k) 
distribution  of  Definition  2.5  and  Srinivasan  (1970),  with  ft"  * 

(Exp(B) :  B  >  0}. 

[Note:  In  Section  11,  it  will  be  seen  that  the  T,  used  in  the  two 
decision  rules  above,  is  the  MLE  of  B.] 
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Similar  decision  rules  exist  for  Type  II  Censoring  and  for  other 
goodness-of-fit  statistics.  These  rules  will  not  be  given  here.  The 
given  rules  are  illustrated  by  numerical  examples  in  the  appendix. 

The  other  detection  problem  with  two  nuisance  parameters  to  be 
considered  here  is  the  two-sample  problem. 

7.  Signal  Detection:  A  Two-Sample  Problem. 

Here,  one  is  comparing  two  sets  of  data,  and  the  PN-situation  is 
the  equality  of  the  two  Weibull-type  laws. 

PN2:  Lj  *  L2  vs.  N  +  S2:  Lj  /  I2 

The  decision  rule  here  will  be  given  for  the  case  of  Type  II  Censoring 

only.  However,  the  corresponding  rule  for  Type  I  Censoring  is  very 

similar  and  can  be  easily  derived. 

Consider  data  Z  =  [W, ,  . . . ,  W  ,  V, ,  ...»  V  ]  and  Y's  and  V's 
<v,  l  ml  n 

as  in  Theorem  5.3. 

Decision  Rule  7.1.  Decide  N  +  S-  iff  —  <  f  ’  or  >  f",  where 

- £  U 

f*  and  f"  are  the  appropriate  percentiles  of  an  Ff2m,  2n)  distribution. 

The  latter  decision  rule  is  based  on  the  minimal  sufficient 
statistics  of  the  Y-sample  and  U-sample  and  is  presumed  to  be  optimal. 
Therefore  no  other  decision  rule  will  be  considered  for  this  case. 

The  next  detection  problem  involves  no  nuisance  parameters. 
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8.  Signal  Detection:  Goodness -of- fit. 

One  has 

PN3:  L  e  fi(WPP)  and  (a, 6)  =  (o0>  8Q)  vs 
N  +  S3:  Not  PN3 


Decision  Rule  8.1.  Decide  N  +  S3  iff  N(T*)  £  Cj  or  C2 


or  if  Cj  <  k  =  N(T*)  <  C2  and 


.  -  R.Z 

D(k)  =  sup  |±-  l  £  (z  -  Y.)  -  [1  -  e  ]|  >  d(a,k2) 


*k+1, 

where  =  -In  — j*  and  and  C2  are  appropriate  percentiles 

B0 

of  a  Poisson  distribution  with  mean  V  =  a^fT*) 

The  decision  rule  above  involves  both  a  "size"  statistics  N(T*) 

and  a  shape  statistic  D(k).  The  PN- distribution  of  N(T*),  the  total 

B0 

number  of  waiting  times ,  is  Poisson  with  parameter  v  =  aQ  [T*]  and 
does  not  involve  the  distribution  of  waiting  times*within  the  interval 
fO,  T*].  On  the  other  hand,  the  statistic  D(k)  depends  directly 
on  the  distribution  of  the  k  waiting  times ,  and  is  hence  a  "shape" 
statistic. 

These  statistics  are  used  separately  in  the  one-parameter  problems 


discussed  below. 
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9.  Signal  Detection:  One-parameter  Problems. 

There  are  two  basic  one-parameter  detection  problems  for  Q(WPP). 

1 

PN^:  0  =  8q  vs.  N  +  S4:  6  ^  8Q;  and 

PNj. :  a  =  a0  (assuming  6,  known)  vs.  N  +  S(.:  a  aQ. 

For  the  first  problem  above,  one  utilizes  the  statistic  D(k) 
of  Decision  Rule  8.1.  For  the  second  problem  above,  one  employs  the 
statistic  N(T*)  of  the  same  decision  rule. 

It  is  interesting  to  note  that  one  need  not  know  a,  in  order 
to  examine  questions  about  8;  but  the  reverse  is  not  true. 

One  now  turns  to  estimation  problems,  where  the  emphasis  will  be 
on  MLE's,  unbiased  point  estimates  and  confidence  intervals. 


10.  Estimation  when  the  Shape  Parameter,  8,  is  known. 


Theorem  10.1.  The  likelihood  functions  for  the  four  sampling  plans 
are  as  follows.  (See  Section  3  for  a  description  of  the  plans  and 
associated  data.) 

(1)  Plan  I:  L(k;  Wj,  ....  w^)  =  (a8)k(ir  w^)®'1  exp  {-a(T*)®>. 


(2)  Plan  II:  L(w v  ....  wk)  =  (a6)k(wk)e_1(Vw.)e'1  exp  {-a(wk)6}. 

1  k 

k  lYi 

(3)  Plan  III:  L^,  ....  yk)  -  Ory^l)-1  exp  {-kot®}  •  (at®)1  , 
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where  Y_.  =  N(t^)  -  N(t^_j);  I  Yj  =  N(tk)  and  tk  * 

k 

IXj  k 

(4)  Plan  IV:  L(yx>  ....  yk)  =  C^~D 1  (^y)  .  where 


Y.  =  N(w?)  -  N(w?  ,). 
3  3  3- 1 


Now  let  a  be  the  MLE  for  the  plan  under  consideration;  and 
a  be  the  associated  unbiased  estimate  of  g. 

These  statistics  and  their  distributions  are  given  below. 

Theorem  10.2.  The  MLE's,  unbiased  estimates  and  sampling  distributions 
for  the  four  plans  are  as  given  in  Table  10.1  below. 

TABLE  10.1.  Estimating  a  when  8  is  known. 

(MLE's;  unbiased  Estimates  and  Distributions) 


Sampling 

a,  the  associated 

Relevant 

Plan 

a,  the  MLE 

unbiased  estimate 

Distribution 

I 

(Type  I 
Censoring) 

N(T*) 

er*r 

A 

a 

(T*)6£  M>(v). 

V  g 

with  v  =  a(T*) 

II 

k 

C  k  )a 

(Type  II 

(wk)8 

”  *  x2k 

Censoring) 

u 

III 

N(W*) 

ka  'v  NB(k,  p) 

(Same-Shape 

Is 

fc 

a 

where  d  *  — — 

Sampling) 

a+i 

IV 

»<v 

.  6a  S'' 

kt.a  *  t,a 

(Same- Dis¬ 

<./ ; 
N(txk1/8) 

1  k 

tance  Samp¬ 
ling) 

a 

'V'  PQ(v),  with 
v  ■  at® 

* 

v  -  atk 
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Confidence  interval  for  a  can  now  be  constructed  using  the 
relevant  distributions  of  Table  10.1 

Definition  10.1  Let  X  'v  X^;  N*  'v  PQ ( v) ;  and  U*  ^  W  -  8(r,p). 

Let  (b^m),  b2(m)],  [Cj(N*),  C2(N*)],  and  fd^r.p),  d2(r,p)]  be 

defined  by  1  -  y  -  Ptb^On)  <_  X  £  b2(m)}  =  P{C^(N*)  <_  V  £  C2(N*)}  = 
PfdjCU'.r)  <  p  <  d2(U*,r)} 

[Note;  For  the  c-values,  see  Rao  et.al . (1966) ;  and  for  the  d-values, 
see  Clemans  (1959)]. 

The  confidence  intervals  of  level  (1  -  Y)  for  a  can  be  constructed 
as  in  the  theorem  below. 

Theorem  10.3.  Die  following  are  respective  confidence  intervals 
of  levels  (1  -  y)  for  a 

(1)  PLAN  I: 

(Type  I  Censoring) 

(2)  PLAN  II: 

(Type  II  Censoring) 

(3)  PLAN  III: 

(Same-Shape) 

(4)  PLAN  IV: 

(Same -distance) 


with  the  indicated  sampling  plans. 


C,(N(T*))  C  (N(T*)) 

1  <  a  <  - 


(T*T 


(T*)6 


ab  (2k)  ab  (2k) 

<  a  < 


2k 


2k 


dj(N(w*),k)  d2(N(w*),k) 

1-dj  (N(w£)  ,kT  <  a  <  :i-d2(N(w*>,k) 


C,(N(tv))  C  (N(tk)) 

1  *  <  a  <  — 


<v 


6 


6 
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One  notes  that  N(T*)  and  N(t^)  are  discrete  random  variables, 
and,  hence,  only  a  discrete  set  of  confidence  levels  is  attainable 
for  Plans  I  and  IV. 

Numerical  examples  illustrating  these  techniques  are  presented  in 
the  appendix. 

When  8  is  not  known,  the  estimation  problems  are  quite  different. 
In  fact.  Plans  III  and  IV  cannot  be  implemented,  and,  hence,  one 
considers  only  Type  I  and  Type  II  censoring. 

11.  Estimation  when  3  is  unknown 

For  Type  I  censoring,  one  can  prove 

Theorem  11.1.  (1)  The  MLE's  of  a  and  8  satisfy  the  equations 

_ N(T*) _ 

k 

N(T*)ln  T*  -  l  InW. 

1  3 


Recalling  the  definitions  of  the  previous  section,  one  has 

Theorem  11.2. 

b  (2k)8  b  (2k)g 

P  {2NCFT  <  8  <  2N(T*)  *  1  "  y 

This  confidence  interval  does  not  depend  on  a.  Further,  there 
is  no  simple  tractable  confidence  interval  available  for  o. 


n  -  NCT*1 

’  (1*)B 


and 


(2) 


2BN(T*) 

A 

6 


%  X 


2k 
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For  Type  II  censoring,  it  can  be  proved  that 


Theorem  11.3.  (1)  The  MLE's  of  a  and 

A  V  A  V 

«  -  -—r  and  e - — —  ; 

(Wk)8  (k-l)lnWk  -  l  W. 


(2) 


'u  X 


2 

2 (k— 1) 


B  satisfy  the  equations 
and 


The  resulting  confidence  interval  for  B  is  given  by 

b  (2(k-l))B  b  (2 (k-l))B 

Theorem  11.4.  P  {  - -  <  B  <  - ~ -  }  =  1  -  Y. 


Again,  no  tractable  confidence  interval  for  a  is  available  in 
this  case. 

The  final  section  discusses  conclusions  and  open  problems. 


12.  Concluding  Remarks;  Open  Problems 

(A)  Applications ■  It  is  clear  that  fi(WPP)  contains  a  wide 
variety  of  NHPP's,  and,  therefore  can  be  used  to  model  a  variety  of 
types  of  data.  From  the  problems  posed  and  solved,  it  is  also  clear 
that  ft(WPP)  is  a  rich  source  of  problems  -  both  theoretical  and 
applied. 

(B)  Other  Sampling  Plans.  Only  two  sampling  plans  were  considered 
when  B  is  unknown.  It  is  desirable  to  develop  other  reasonable  nlans  , 
since  in  applications  B  will  rarely  be  known.  Regularly  spaced 
sampling  yields  data  of  the  form  N  ■  [N(A) ,  N(2A),  ...,  N(kA)].  But 

'V 

the  MLE's,  etc.,  are  totally  intractable. 
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(C)  Optimal  Procedures .  For  each  of  the  inference  problems 
considered,  one  or  more  reasonable  procedures  is  given.  Are  these 
procedures  optimal?  Should  goodness -of-fit  statistics  other  than 
the  K-S  statistics  be  used.  For  example,  is  the  Cramer-von-Mises 
statistic  better  here? 

(D)  Confidence  Intervals  for  a.  When  8  is  unknown,  the  authors 
were  unable  to  construct  reasonable  confidence  intervals  for  a. 

Would  another  type  of  sampling  plan  yield  better  results  in  this 
direction? 

(E)  Joint  Confidence  Regions .  Does  the  distribution  theory  pre¬ 
sented  allow  one  to  construct  reasonable  confidence  regions  for  the 
pair  (a, 8)?  This,  for  example,  is  readily  done  for  certain  Gaussian 
processes . 

(F)  Welbull  Fits.  In  the  body  of  the  text  a  goodness-of-fit 
test  is  given  for  a  WPP.  However,  in  seeking  to  model  engineering  or 
biomedical  data,  etc. ,  with  a  WPP,  it  would  be  interesting  to  develop 
a  set  of  axioms  which  the  underlying  physical  forces  should  satisfy  in 
order  to  generate  a  WPP.  Such  axiom  systems  are  common  for  HPP's, 
for  example.  Again,  one  asks  about  the  feasibility  of  using  goodness- 
of-fit  tests  other  than  the  K-S. 
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APPENDIX:  NUMERICAL  EXAMPLES 

In  the  appendix  one  is  concerned  with  three  sets  of  data  in 
illustrating  the  methodology  developed.  Table  A.l.  contains  F.KG- 
data  for  25  cycles  of  a  "normal"  heart.  Table  A.2.  contains  simulated 
WPP  data  with  several  sets  of  parameter  values.  In  the  n»  ^rical 
examples.  Data  Sets  2  and  3  are  used  (as  well  as  the  EKG-data) . 

(I)  Numerical  Examples  for  Section  6.  There  are  three  decision 
rules  in  this  section,  and  they  are  illustrated  using  the  data  of 
Table  A.l.  Modifications  of  the  data  to  be  used  with  the  various 
decision  rules  of  Section  6  are  given  in  Tables  A. 2,  A. 3,  and  A. 5. 


i 

i 

\ 


TABLE  A. 1 .  EKG-Data 
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MAGNITUDE  TIMES 


p 

Q 

R 

T 

P 

Q 

—  —  ■ 

R 

T 

Cycle 

Length 

.0927 

.0708 

1.190 

.  4785 

148 

428 

B 

1014 

.1005 

.0917 

1.203 

.4344 

1,183 

1,440 

1,484 

1,758 

1014 

.1121 

.0982 

1.212 

.4316 

1,897 

BBSS 

2,498 

HN 

1040 

.1238 

.0608 

1.297 

.  5652 

3,452 

3,521 

mm 

WB 

1041 

.0819 

.0206 

1.261 

.  5060 

4,108 

m 

1056 

.0769 

.0764 

1.284 

BUS 

5,346 

suss 

5,655 

IBB 

1078 

.0815 

.0507 

1.228 

.  4  7  G  0 

6,346 

6,689 

RS 

:  1073 

.1263 

.1113 

1.166 

.  4240 

7,596 

7,806 

8,085 

1076 

.1009 

.0817 

1.274 

Kffl 

8,645 

8,840 

8,882 

9,154 

1089 

.0839 

.0601 

1.239 

.4663 

9,694 

9,935 

m 

10,250 

1044 

.1126 

.0907 

1.164 

.4462 

10,727 

11,015 

1012 

.0939 

.0894 

1.189 

m 

iws 

11,985 

12,027 

1028 

.0929 

.0648 

1.278 

Rum 

12,876 

13,055 

13,331 

1038 

.1010 

.0555 

1.213 

.4742 

IBHBI 

m 

mm 

BBI 

1062 

“sra — 

.0560 

.0374 

1.235 

.4874 

mm 

HUB 

15,153 

HB 

1063 

.0832 

.0899 

1.  294 

.4268 

16,179 

16,216 

16,495 

1083 

.1278 

.0939 

1.  294 

.4465 

BBSS 

17,265 

17,299 

17,587 

1063 

.1044 

.0585 

1.255 

.  5001 

18,255 

18,362 

1018 

.1184 

.0702 

1.250 

.4491 

EBBS 

EBBS 

19,380 

19,655 

1013 

.1142 

.1141 

1.266 

.  4105 

20,203 

EB9S' 

947 

.1338 

.1140 

1.228 

.4052 

^^9 

hssu 

MB 

891 

.1331 

.1050 

1.  230 

.4405 

mm 

911 

ITS - 

.0849 

.0584 

1.127 

.  4265 

mill 

901 

.0916 

.0807 

1.156 

.4038 

sn 

gfgl 

mm 

947 

.1442 

.1323 

1.260 

i 

i 

Bill 

ill 

— - 
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TABLE  A. 2.  Illustration  of  Decision  Rule  6.1. 


X. 

J 


W.  Y. 

J  J 


Z 


r 


E 

r 


E(j)- 


JdL 


k-1 


k-1 

E(j) 


1 

1014 

1014 

0.0200 

0.46 

0.0035 

0.0035 

0.0441 

2 

1014 

2028 

0.0590 

1 . 758 

0.0133 

- 

0.0776 

3 

1040 

3068 

0.1000 

3.858 

0.0291 

- 

0.1073 

4 

1041 

4109 

0.1419 

6.696 

0.0506 

- 

0.1312 

5 

1056 

5165 

0.1883 

10.2737 

0.0776 

- 

0.1497 

6 

1078 

6243 

0.2406 

14.6045 

0.1103 

- 

0.1624 

7 

1073 

7316 

0.2960 

19.6365 

0.1484 

- 

0.1698 

8 

1076 

8392 

0.3573 

25.3533 

0.1915 

- 

0.1721 

9 

1089 

9481 

0.4239 

31.7118 

0.2396 

- 

0.1695 

10 

1044 

10525 

0.4938 

38.625 

0.2918 

- 

0.1627 

11 

1012 

11537 

0.5690 

46.022 

0.3477 

- 

0.1523 

12 

1028 

12565 

0.6484 

53.8028 

0.4065 

• 

0.1390 

13 

1038 

13603 

0.7337 

61.8735 

0.4675 

- 

0.1234 

14 

1062 

14665 

0.8255 

70.1285 

0.5298 

- 

0.1066 

15 

1063 

15728 

0.9300 

78.4985 

0.5931 

- 

0.0887 

16 

1083 

16811 

1.0520 

86.9145 

0.6566 

- 

0.0707 

17 

1063 

17874 

1.1892 

95.2389 

0.7195 

- 

0.0532 

18 

1018 

18892 

1.3478 

103.3257 

0.7806 

0.0079 

0.0376 

19 

1013 

19905 

1.5374 

111.0127 

0.8387 

0.0205 

0.0249 

20 

947 

20852 

1.7661 

118.0771 

0.8921 

0.0285 

0.0170 

21 

891 

21743 

2.0583 

124.252 

0.9387 

0.0296 

0.0158 

22 

911 

22654 

2.4723 

129.1966 

0.9761 

0.0216 

0.0239 

23 

901 

23555 

3.1654 

132.362 

24 

947 

24502 

24  Data  points  T*  =  24030,  k  =  23 

Let  {Xj}  be  the  interarrivals  and  {W^ }  be  the  waiting  times. 


Also  let  Y_.  =  -1 
2r 

and  E  =  ■=—  . 


=  -in  (-*£!),  z  =  l  (k+l-j)Y. 
3  T  r  j-1  3 


k  W 

-  I  a  ln(=£) 
m=k+l-r 


The  critical  value  (on  interpolating)  is  0.254.  is  computed 

.  Therefore,  one  decides  PN. 


to  be  0.1721 
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TABLE  A. 3.  Illustration  of  Decision  Rule  6.2. 


Y. 

J 

*  =l-exp  {- 

Y 

*  .  hi 

23 

23  '  * 

1 

0.0200 

0.0221 

0.0221 

0.0214 

2 

0.0590 

0.0637 

- 

0.0280 

3 

0.1000 

0.1056 

0.0186 

0.0248 

4 

0.1419 

0.1464 

0.0160 

0.0275 

5 

0.1883 

0.1895 

0.0156 

0.0279 

6 

0.2406 

0.2354 

0.0180 

0.0255 

7 

0 . 2960 

0.2812 

0.0203 

0.0231 

8 

0.3573 

0.3288 

0.0245 

0.0190 

9 

0.4239 

0.3768 

0.0290 

0.0145 

10 

0.4938 

0.4236 

0.0323 

0.0112 

11 

0.5690 

0.4700 

0.0155 

0.0083 

12 

0.6484 

0.5149 

0.0366 

0.0068 

13 

0.7337 

0.5589 

0.0372 

0.0063 

14 

0.8255 

0.6019 

0.0367 

0.0068 

15 

0.9300 

0.6457 

0.0370 

0.0065 

16 

1.0520 

0.6908 

0.0386 

0.0049 

17 

1.1892 

0.7347 

0.0390 

0.0044 

18 

1.3478 

0.7777 

0.0386 

0.0049 

19 

1.5374 

0.8201 

0.0375 

0.0060 

20 

1.7661 

0.8606 

0.0345 

0.0090 

21 

2.0583 

0.8994 

0.0298 

0.0136 

22 

2.4723 

0.9366 

0.0236 

0.0229 

23 

3.1654 

0.9707 

0.0142 

0.0293 

For  this  data 

Y  =  0.89634  and  DR  =  0.0390. 

Interpolation  in  the  Lilliefors  (1969)  table  yields  a  critical  value 
of  0.22.  The  decision  is  therefore:  PN. 


29  - 


TABLE  A. 4.  Illustration  of  Decision  Rule  6.3. 


Y. 

1 

vr-\  k-1 

*  I Dl) 

kY 

"  23~ 

23  ‘  * 

1 

0.0200 

0.0211 

0.0211 

0.0224 

2 

0.0590 

0.0611 

0.0176 

0.0259 

3 

0.1000 

0.1015 

0.0145 

0.0279 

4 

0.1419 

0.1410 

0.0106 

0.0329 

5 

0.1883 

0.1828 

0.0089 

0.0346 

6 

0.2406 

0.2276 

0.0102 

0.0333 

7 

0 . 2960 

0.2725 

0.0116 

0.0318 

8 

0.3573 

0.3193 

0.0150 

0.0285 

9 

0.4239 

0.3669 

0.0191 

0.0244 

10 

0.4938 

0.4134 

0.0221 

0.0214 

11 

0.5690 

0.4598 

0.0250 

0.0185 

12 

0.6484 

0.5049 

0.0266 

0.0168 

13 

0.7337 

0.5494 

0.0277 

0.0158 

14 

0.8255 

0.5930 

0.0278 

0.0157 

15 

0.9300 

0.6378 

0.0291 

0.0144 

16 

1.0520 

0.6841 

0.0319 

0.0116 

17 

1.1892 

0.7294 

0.0337 

0.0097 

18 

1.3478 

0.7741 

0.0350 

0.0085 

19 

1.5374 

0.8182 

0.0356 

0.0079 

20 

1.7661 

0.8606 

0.0345 

0.0090 

21 

2.0583 

0.9011 

0.0315 

0.0119 

22 

2.4723 

0.9398 

0.0268 

0.0167 

23 

3.1654 

0.9744 

0.0179 

0.0256 

One  computes  Y  =  0.89634 

=  0.0356  and  the  internolated  Srinivasan  value 


is  0.15.  Therefore,  one  decides:  PN 
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TABLE  A. 5.  Simulated  Waiting  Times  with  y(t)  =  at 


Data  Set  1 

Date  Set  2 

Date  Set  3 

Data  Set  4 

a  =  4.0 

a  =  4.0 

a  =  2.5 

a  =  1.0 

$  =  0.5 

e_=i_,o _ 

6  =  2.0 

8-2.0 

1 

0.179ii 

0.261*0  * 

0.6621 

1.5370 

2 

0.2118 

0.3535 

0.8156 

1.81*15 

3 

3.511*3 

0.8270 

1.1033 

2.2670 

1* 

3.5301 

1.1*313 

1.11*25 

2.1*329 

5 

1>.8938 

1.7190 

1.3875 

2.6718 

6 

5.1777 

2.3331* 

1.1*569 

2.9192 

7 

5.2139 

2.7005 

1.6535 

2.9655 

8 

5.U513 

2.9981 

1.6551* 

3.3051 

9 

5.5830 

3.031*9 

1.7187 

3.3359 

10 

5.8101 

3.1*578 

1.8768 

3.6561 

11 

6.6792 

3.8193 

2.121*2 

3.6623 

12 

7.1621* 

lt.1066 

2.1813 

3.81*65 

13 

7.1868 

1*.  281*1 

2.2706 

3.8716 

1U 

8.3186 

1*.2980 

2.3851 

3.9012 

15 

10.1831* 

lt.3277 

2.1*161 

1*.  1*260 

16 

17.71*97 

1*.  1*119 

2.5016 

1*.  1*830 

17 

20.938!* 

1*.9939 

2.5161 

U.6552 

16 

37.81i07 

5.7613 

2.5183 

1*.  71*08 

19 

1*2.9353 

6.1637 

2.5231 

1*.  8U1I 

20 

1*7.5891 

6J*9i*8 

2.571*5 

1**8968 

21 

1*7.6528 

6.6561 

2.7927 

it. 9901* 

22 

1*8.2519 

7.171*8 

2.9971 

5.0636 

23 

1*8.8106 

7.21*00 

3.0902 

5.0779 

2l* 

1*9.7105 

7.3192 

3.131*6 

5.0933 

25 

53.5510 

7.1*199 

3.1652 

5.21*79 

26 

51*.  3202 

7.7230 

3.1985 

5.1*329 

27 

55.3261 

7.7686 

3.2867 

5.1*1*93 

28 

56.291*9 

7.9861* 

3.3588 

5.1*559 

29 

67.6298 

8.3356 

3.1*81*9 

5.5021 

30 

72.0625 

9.1*819 

3.5187 

5.6157 

31 

72.1961 

9.7305 

3.51*76 

5.6721 

32 

75.2956 

10.1198 

3.5925 

5.71*19 

33 

78.5233 

10.3810 

3.8160 

5.9278 

31* 

79.6223 

10.5817 

3.81*90 

5.9652 

35 

81.8531 

10.811*2 

3.9115 

6.01*53 

36 

89.81*1*9 

10.9833 

3.9896 

6.1611 

37 

95.1*882 

11.1056 

U.0751* 

6.2578 

38 

100.9558 

11.3623 

U.1001* 

6.3962 

39 

108.311*0 

11.6609 

ll.ll85 

6.3982 

i*o 

115.3873 

11.7117 

L.1229 

6.1*602 

la. 

122.6267 

12.1708 

1*.2126 

6.5531 

1*2 

139.9208 

12.9102 

It.  291*9 

6.5856 

1*3 

11*2.7151 

12.9625 

lt.2966 

6.5952 

1«1* 

11*6.1*252 

13.6938 

1».3525 

6.6012 

1*5 

11*6.61*17 

Ui.5538 

1*.  1*296 

6.6321 

1*6 

11*8.981*8 

15.2888 

1*.  1*329 

6.7129 

1*7 

205.9331 

15.3991* 

1*.  1*1*72 

6.7171 

1*8 

211.0332 

15.1*766 

1*.  51*33 

6.9320 

1*9 

217.9563 

15.5336 

It.  5572 

7.0581 

50 

239tJ*03U 

15.5897 

It.  5681 

7.181*9 

51 

250.9719 

16.3861* 

1*. 5883 

7.2028 
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TABLE  A. 6.  Illustration  of  Decision  Rule  7.1. 

[Here  one  uses  simulated  Data  Set  2  with  a  =  4.0,  6  »  1.0] 

The  first  24  points  constitute  the  first  sample;  and  the  next  20  points 
constitute  the  second  sample. 


1st  Sample  2nd  Sample 


W. 

Y. 

V. 

U. 

3 

3 

3 

3 

1 

0.2640 

0.0137 

1 

0.3031 

0.1429 

2 

0.3535 

0.0245 

2 

0.3487 

0.2524 

3 

0.8270 

0.0336 

3 

0.5665 

0.2619 

4 

0.4313 

0.1086 

4 

0.9157 

0.4065 

S 

0.1790 

0.1332 

5 

2.062 

0.5082 

6 

2.3334 

0.1855 

6 

2.3106 

0.5201 

7 

2.7005 

0.2530 

7 

2.6999 

0.5931 

8 

2.9981 

0.3959 

8 

2.9611 

0.6604 

9 

3.0349 

0.5199 

9 

3.1618 

0.6941 

10 

3.4578 

0.5391 

10 

3.3943 

0.7428 

11 

3.8193 

0.5460 

11 

3.5634 

0.8137 

12 

4.1066 

0.5493 

12 

3.6857 

0.8793 

13 

4.2841 

0.5916 

13 

3.9424 

0.9716 

14 

4.2980 

0.6641 

14 

4.241 

1.1204 

15 

4.3277 

0.7635 

15 

4.2918 

1.2412 

16 

4.4119 

0.8940 

16 

4.7509 

2.0529 

17 

4.9939 

0.9062 

17 

5.4903 

2.5331 

18 

5.7613 

1.0107 

18 

5.5426 

3.0184 

19 

6.1637 

1.1568 

19 

6.1839 

3.1586 

20 

6.4948 

1.4624 

20 

7.1339 

21 

6.6561 

1.6456 

22 

7.1748 

2.1941 

23 

7.2400 

3.0440 

24 

7.3192 

3.3360 

25 

7.4199 

One  computes  V*  *  7.1339,  W*  =  7.4199,  Y  =  0.8738,  U  =  1.0831 

Y 

and  —  =  0-8068.  The  interpolated  F-values  for  48  and  38  degrees 
U 

of  freedom  are  0.5473  and  1.8676.  Therefore,  one  decides:  PN. 


Numerical  Examples  for  Sections  8  and  9. 


The  Decision  Rule  8.1  is  a  combination  of  two  decision  rules. 


The  application  of  these  two  rules  is  given  below  in  the  context  of 
Section  9  and  the  EKG-Data  of  Table  A.l. 


One  first  considers 


PN.:  6  ■  B„  vs.  N  ♦  S.:  B  i  Bn 

(S0-2)  0  - 0 

0  i  k  -v 

The  statistics:  D(k)  «  sup  l  e  (z  -  Y.)  -  [1  -  e  ]| 

z  1  J 

W 

where  Y.  =»  -In  [-—^1  ,  (T*  =  11031)  . 

1  T 
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TABLE  A. 7.  Illustration  relative  to  PN,  of  Section  9. 

-  4 


X. 

J 

W. 

J 

Y. 

J 

*  ,  "0oz 

*  =  1-e 

*_  tl 

10 

i_ . 
10 

1 

1014 

1014 

0.0470 

0.0897 

0.0897 

0.0103 

2 

1014 

2028 

0.1514 

0.2613 

0.1613 

- 

3 

1040 

3068 

0.2734 

0.4212 

0.2212 

- 

4 

1041 

4109 

0.4106 

0.5601 

0.2601 

- 

5 

1056 

5165 

0.5692 

0.6797 

0.2797 

- 

6 

1078 

6243 

0.7588 

0.7808 

0.2808 

- 

7 

1073 

7316 

0.9875 

0.8612 

0.2612 

- 

8 

1076 

8392 

1.2797 

0.9226 

0.2226 

- 

9 

1089 

9481 

1.6937 

0.9662 

0.1662 

- 

10 

1044 

10525 

2.3868 

0.9916 

0.0916 

0.0084 

One  computes  D(k)  =  0.2808,  and  finds  the  critical  value  to  be  ■  0.369. 
Therefore,  one  decides  PN. 

For  the  next  illustration,  consider 

PNg :  a  *  .072  vs.  N  ♦  Sg:  a  t  .072 

(Assuming  3  known  to  be  0.4) 

Decision  Rule.  Decide  N  +  S  iff 

N(  *)  >  8  or  <0 

Here  the  FAR  %  .05 
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We  know  that  N(T*)  %  PQ(V)  where  V  =  ct(T*)6  =  ( .072)  (11031) A  = 

2.98. 

From  the  data,  N(T*)  =  10.  Therefore,  one  decides  N  +  S,.. 


(IV)  Numerical  Examples  for  Section  10. 

The  data  is  that  of  Data  Set  3  of  Table  A. 6. 

(A)  The  point  estimates  from  Theorem  10.2  are  as  follows: 


PLAN  I: 


=  _  '  _  N(T*) 
a  =  a  =  - - 

(T*) 


50 

(4. 58)  2 


2.384 


(Here,  take  T*  =  4.58  and  0=2.  From  the  data  one  sees  N(T*)  *  50). 


PLAN  II:  (a)  a 


k  =  51. 


51 


(Wk)p  (4.5833) 


Y~  =  2.4278,  when  one  chooses 


'u  ir  i  en 

(bj  a  =  — a  =  (fy)  (2 >4278)  =  2.3802 


PLAN  III: 


a  =  a 


N(W*) 

~T~ 


47 

16 


2.9375. 


Here,  one  chooses  k  =  16;  and  generated  the  first  16  waiting  times 

2 

{W*}  of  a  WPP  with  y(t)  =  t  .  It  turns  out  that  for  this  second 
n 

process,  which  is  independent  of  the  process  of  interest,  W*^  =  4.4830 
and  N(W*6)  *  47. 


I 


i 


! 


I 


i 

i 


i 
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*  %  N(t,  ) 

PLAN  IV:  a  =  a  -  - 5—  =  - -  =  2.384  when  one  chooses 

(tk)p  (4.58j 

tk  =  4.58  and  k  =  16.  This  means  that  =  (tk)(k)-1^  =  1.145 
and  t.  =  tjCj)176  =  (1.145)  /j  for  j  =  1,  2 .  16. 

(B)  The  Confidence  Intervals  based  on  Theorem  10.3. 

PLAN  I;  (Type  I  Censoring).  Here  0  =  2.0  is  known;  T*  =  4.58,  k  =  50 
from  Data  Set  3  and  one  seeks  a  95%  confidence  interval  for  a. 

One  has  N(T*)  n,  PQ(v),  where  v  =  a(T*)B  *  a(4.58)2.  From  the  table 
on  p.  46  of  Rao  et. al. (1966) ,  one  finds 

Cj(N(T*))  C2(N(T*)) 

[ - - -  ,  - g - ]  leads  to 

(T*)p  (T*)p 

— 57,67  =  1.7958  <  a  <  3.0963  =  64,95 
(4.58j  (4. 58) 2 

PLAN  II:  (Type  II  Censoring).  Here  k  «  51,  Wk  =  4.5833  and  0  *  2.0 

2 

is  known;  and  one  needs  the  X  -percentiles  for  102  degrees  of  freedom 

Sbx(2k)  Sb2(2k) 

[ — 2lc -  »  — 2k - J  leads  to 


(2.4278) (76. 25 
(2) (51) 


1.8149  <  a  <  3.1547 


(2.4278) (132. 54 
2(51) 
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PLAN  III:  (Same-Shape)  Here  one  uses  Data  Set  3,  and  an  independent 

8  2 

WPP  with  mean  function  p(t)  =  tP  =  t  .  For  k  =  16,  W*  =  4.4830 

and  N(W*)  =  47.  Since  N(W£)  has  a  negative  binomial  distribution, 

one  employs  the  methodology  and  graphs  of  Clemans  (1959),  to 

construct  a  90%  confidence  interval  for  the  mean,  m*,  of  this  negative 

.  N(wk*)  47 

binomial  distribution,  based  on  m*  =  — ^ —  =  jg  =  2.9375.  The  90% 
band  on  m*  from  Clemans  graphs  is  2.1  <  m*  <  5.2.  But  in  terms 

a 

of  o,  m*  =  ■  °y* -  »  a,  and  hence  2.1  <  a  <  5.2. 

a+T 


PLAN  IV  (Same- Distance)  Here  one  again  uses  Data  Set  3;  assumes  $  -  2 
is  known;  and  chooses  tj  =  1.07.  Then  tj  =  (tj)(j)*^  =  1.07  /f,  for 
j  -  1,  2,  ....  16 


Cj(N(tk)) 


=  1.7958  <  a  <  3.0963  = 


C2(N(tk)) 

~ivr 


One  now  considers  the  cases  when  8  is  unknown. 


(V)  Illustrations  of  the  Methodology  of  Section  11. 


PLAN  I.  (Type  I  Censoring)  Here  one  uses  Data  Set  I  of  Table  A. 6,  and 
T*  *  245.  (The  formulae  are  from  Theorems  11. 1  and  11.2). 


(a)  Point  Estimates:  The  MLE's  are 

50 


8 


N(T*) 


and  J-® 
(T*)P 


(245) 


0.4894 


3.386. 


N(T*)lnT*  -  l  InW 
1  3 


4894. 
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(b)  95%  Confidence  Interval  for  g. 


bj (2k) § 

2N(T*) 


=  0.3648  <  6  <  0.6325  = 


b2(2k)B 

2N(T*)'" 


2 

where  the  X  -percentiles  (b's)  are  approximately  74.55  and  129.24 
for  100  degrees  of  freedom. 


PLAN  II  (Type  II  Censoring)  (The  formulae  are  from  Theorems  11.3  and 
11.4).  One  uses  Data  Set  I  and  k  =  5,  from  which  W51  =  250.9719. 

(a)  Point  Estimates.  One  easily  calculates 


B  = 


k-1 

(k-l)ln  W.  -  l  W. 


0.4934  and 


a  = 


51 


(250.9719) 


0.4934 


=  3.273. 


(b)  95%  Confidence  Interval  for  B. 

Tbe  b's  are  as  in  the  illustration  of  PLAN  I,  i.e.,  bj(lOO) 
74.55  and  b^lOO)  =  129.24,  interpolating  from  a  X^-table.  Hence 

b,  (100)6  b  (100)6 

— - -  =  0.3606  <  B  <  0.6252  =  — - - 


This  completes  the  numerical  examples. 


